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ABSTRACT 

The current dynamo paradigm for the Sun and sun-like stars places the generation site for strong 
toroidal magnetic structures deep in the solar interior. Sunspots and star-spots on sun-like stars are 
believed to arise when sections of these magnetic structures become buoyantly unstable and rise from 
the deep interior to the photosphere. Here we present the first 3-D global magnetohydrodynamic 
(MHD) simulation in which turbulent convection, stratification, and rotation combine to yield a 
dynamo that self-consistently generates buoyant magnetic loops. We simulate stellar convection and 
dynamo action in a spherical shell with solar stratification, but rotating three times faster than the 
current solar rate. Strong wreaths of toroidal magnetic field are realized by dynamo action in the 
convection zone. By turning to a dynamic Smagorinsky model for subgrid-scale turbulence, we here 
attain considerably reduced diffusion in our simulation. This permits the regions of strongest magnetic 
field in these wreaths to rise toward the top of the convection zone via a combination of magnetic 
buoyancy instabilities and advection by convective giant cells. Such a global simulation yielding 
buoyant loops represents a significant step forward in combining numerical models of dynamo action 
and flux emergence. 

Subject headings: convection — magnetohydrodynamics (MHD) — stars: interiors — stars: magnetic 
field — Sun: activity 



1. CONVECTION, ROTATION AND MAGNETISM 

The clearest signature of the global solar dynamo is 
the emergence of sunspots at the photosphere. Creating 
these coherent magnetic structures likely requires sev- 
eral dynamical processes operating at various locations 
in the solar interior. A single 3-D numerical simulation 
of solar magnetism that extends from the deep interior 
through the Sun's upper atmosphere, while resolving all 
relevant scales, is intractable with current computational 
resources. This leads to three main classes of simula- 
tions that address elements of solar-like dynam o pro- 
cesses (see reviews iFan 20091 : iCharbonneau 2010h . One 
approach to study how loops may emerge is to insert 
a compact magnetic field structure into a spherical do- 
main and track its buoyant rise (e.g..lCaligari et al 1991 
IFan 20081 Uouve fe Brun 20091 : IWeber et al. 20111 ). An- 
other approach uses local planar models with me- 
chanical forcing to generate large-scale shear that 
drives dyna mo action and c r eates buoyant magnetic 
loops (e.g., iCline et al. 2003t iVasil fe Brummell 20081: 
iGuerrero fc Kapvla, 201l( ). Planar models have also been 
used to study 3-D buoyancy instabilities in a mag- 
netized layer that can lead to rising elements (e.g., 
iKersale et al. 20 071). The third approach uses global 
convective MHD models. These incorporate the rotat- 
ing spherical-shell geometry needed to self-consistently 
generate differential rotation and meri dional circula - 
tion through Reynolds stresses (see review iMiesch 20051 ). 
Such models have captured the formation of mag- 
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netic structures an d cycles in solar (jGhizaru et al. 20101 : 
iRacine et al. 201 il) and rapid l y- rotating sun- like stel- 
lar models (jBrown et al. 20101: [B rown et al. 20lH here- 
inafter B10 and Bll, respectively), yielding differential 
rotation, dynamo action, and large-scale magnetic fields, 
but not buoyant magnetic loops that rise toward the top 
of the convective layer. 

Here we report on a global convective dynamo sim- 
ulation of a sun-like star rotating at three times 
the mean solar angular velocity (3Q0), such as our 
Sun did when it was younger and as do many so- 
lar analogues ([Petit et al. 20 08). This simulation (i) 
attains a differential rotation profile created by the 
interpl ay of convection, r ot ation and stratificatio n 
(e.g., IBrun fc Toomre 20021: IMiesch fc Toomre 20091) . 
(ii) forms global-scale toroidal magnetic structures that 
undergo cycles of magnetic activity and reversals of 
global polarity, and (iii) achieves buoyant magnetic loops 
from the strongest portions of the toroidal structures 
which rise from the base of the convection zone. This 
work extends the work of B10 and Bll in which simu- 
lations of rapidly-rotating suns with moderate levels of 
diffusion were able to accomplish (i) and (ii) . The forma- 
tion of buoyant loops is facilitated in our current work 
by adopting a dynami c Smagorinsky subgrid-scalc model 
(|Germano et al. 19911 ) , which serves to minimize the dif- 
fusion of well-resolved structures. 

2. SIMULATION PARAMETERS AND PROPERTIES 

We have conducted 3-D MHD simulations of turbu- 
lent convection and dynamo action in a spherical shell 
spanning the bulk of the convection zone from 0.72i?© 
to 0.97i? Q involving a density contrast of 25, and ro- 
tating at 3fi Q (1240 nHz, once every 9.3 days). We 
use the anelastic spherical harmonic (ASH) code (e.g., 
IBrun et al. 20041 ) . The anelastic treatment lets us follow 
the subsonic flows in the deep convection zone. Within 



Nelson et al. 




(b) 




(c) 



V r at 0.96 R Q 



B/,at 0.84 Rq 




1190 nHz 135(TnHz -3 k( 



Fig. 1. — Snapshots of flows and fields in case S3 when a buoyant loop begins to rise at time %. (a) Radial velocity v r in global Mollweidc 
projection (equator is dashed) near the top of the computational domain, showing fast, narrow downflows (dark tones) and broad, slower 
upfiows (light tones), (b) Companion view of toroidal magnetic field B^ at mid-convection zone. Effects of turbulent convection contribute 
to the ragged nature of the wreaths. Several buoyant magnetic loops (see Figure [3j are generated in the negative-polarity wreath segment 
just above the equator and right of image center, (c) Time and zonal average of rotation profile (f2), possessing an equatorial region 
with fast rotation and slower rotation at higher latitudes, (d) Longitudinally-averaged toroidal magnetic field (B^) revealing a prominent 
axisymmetric field component. 



this nearly adiabatically stratified region, we expect that 
magnetic buoyancy instabilities captured by our anelas- 
tic treatment differ from fully compressible treatments 
by no more than a few p ercent in either growth rate or 
scale (IBerkoff et al. 20101) . ASH is a large-eddy simula- 
tion (LES) code that resolves the largest scales of motion 
and uses a subgrid-scale (SGS) model to parameterize 
the effects of unresolved, small-scale turbulence. The 
dynamo simulations of B10 and Bll used a SGS model 
where the turbulent magnetic diffusivity rj t was constant 
on spherical shells and in time, and varied only slowly 
with depth as the inverse square-root of the background 
density. B10 examined a simulation (case D3, at 3^0) 
which exhibited persistent toroidal magnetic structures, 
whereas Bll studied a simulation that achieved cycles of 
magnetic activity and global polarity reversals (case D5, 
at SO© ) . These simulations had r/ t = 2.64 x 10 12 and 
1.88 x 10 12 cm 2 s _1 respectively at mid-convection zone. 

Here we consider a new ASH simulation, case S3, which 
achieves much lower levels of diffusion through the use 
of a dynamic Smagorinksy (DSMAG) SGS model. This 
assumes self-similar behavior in the resolved portion of 
the inertial range of scales in a turbulent flow in order 
to extrapolate the effects of unresolved small-scale mo- 
tions on the resolved scales. The resulting viscosity z/g is 
determined by the properties of the grid and the flows, 
and varies by orders of magnitude in all three spatial 
dimensions and in time. To determine the thermal and 
magnetic diffusion coefficients we assume constant ther- 
mal and magnetic Prandtl numbers. In cases D3, D5 and 
S3 these are set to 0.25 and 0.5 respectively. We reserve 
further discussion of the properties of ASH simulations 
using the DSMAG SGS model for a forthcoming paper. 
In case S3 the DSMAG SGS model allows a simulation 
(with 1024 longitudinal, 512 latitudinal, and 193 radial 
grid points) to achieve a mean magnetic diffusion coeffi- 
cient at mid-convection zone of fj t = 4.8 x 10 10 cm 2 s~ 1 . 
This reduction in diffusion by a factor of about 40 from 
case D3 is critical for the formation and coherent rise of 
buoyant magnetic loops. 

Case S3 exhibits turbulent convective patterns shown 
in Figure Ufa) which are largely vortical at high latitudes 
and aligned with the rotation axis near the equator. The 
convection builds and maintains a strong differential ro- 




FlG. 2. — Field reversals with time, (a) Hemispherical volume- 
averaged toroidal magnetic field [B$\ of progenitor case D3b over 
nearly 6000 days, displaying irregular magnetic activity cycles. 
Case S3 branched from case D3b at time t s (dotted line), (b) [BA 
for case S3 over about 800 days. Case S3 continues the cyclic 
behavior of D3b, but additionally produces buoyant loops. The 
creation of loops which pass 0.90-Rq are indicated by tick marks in 
the lower panel. Detailed information on the buoya nt lo op at time 
tb = t s + 683 days (dotted line) is shown in Figures j 1 1 3 1 



tation that is progradc at the equator and retrograde 
at mid to high latitudes (Figure [IJc)). This organized 
shear drives the creation of toroidal magnetic structures 
at low latitudes in each hemisphere, as demonstrated in 
BIO. Here the increased level of turbulence enhances 
the power in smaller-scale components of the toroidal 
field Btj, (Figure Hfb)) while still retaining a substantial 
zonally- averaged toroidal field (B^) (Figure Did)). 

In addition to creating strong magnetic structures near 
the base of the convective region, case S3 also undergoes 
cycles of magnetic activity and reversals of global mag- 



Buoyant Magnetic Loops in a Global Dynamo Simulation 



t b + 4.2d t b +6.0d 



t b +8.2d 




t b + 10.4 d 
A 

- 


t b + 12.6 d 
A 

r 0> 


t b + 14.6 d 

r ^ 




— D 




t 




r 


■ i 1 1 i 1 1 i 1 1 i 1 1 i i 


~, , i . , , , i , , , , i , , 


lllllllll.ll.il! 



Eq +5 +10 Eq +5 +10 Eq +5 +10 Eq +5 +10 Eq +5 +10 Eq +5 +10 Eq +5 +10 

Latitude (degrees) . . 

(b) (c) 

Loop A~" — ■ 



•4 kG 



■40 kG 



0.92 




2 



0.90 



0.88 



73 

CO 
DC 



0.84- 



0.82- 



0.80 



■ Loop A 
* Loop B 

— Advective Motion 

— Mag. Buoy. Motion 



i ■ ■ ■ i ■ ■ ■ i 




r 



0.72 



2 4 6 8 10 12 14 

Time since t b (days) 



Fig. 3. — Analyzing a rising loop, (a) 2-D cuts in longitude at successive times (tracking in longitude at the local rotation rate of the 
loop) showing toroidal magnetic field over radius and latitude. The rising magnetic loop A is seen in cross-section starting at 0.81-Rq at 
t = tf, and rising to 0.91-R© after roughly 15 days. Proto-loop B is also seen rising starting at 8.6 days, but the top of loop B never rises 
above O.88R0. (b) 3-D visualization of magnetic field lines in the core of a wreath which produces four loops (two shown here, one of 
which is loop A) at % + 14.6 days. Perspective is looking down along the rotation axis toward the equatorial plane. Coloring indicates 
field magnitude. Dashed lines indicate radial position. Dotted line shows the cutting plane used in the left-most panel above, (c) Radial 
location of the top of a buoyant loop as a function of time since tf, , along with movement attributable to magnetic buoyancy (red lines) or 
to advection by convective upflows (blue lines). 



netic polarity similar to those described in case D5 in 
Bll. This is consistent with results from parameter sur- 
veys with ASH simulations, which indicate that decreas- 
ing both v and can yield cyclic b ehavior seen at 5£Iq at 
lower rotation rates ((Brown 201 ID . Because of the large 
computational cost of the DSMAG SGS model, case S3 
was started using a less diffusive descendant of case D3 
in BIO, which we labe l case D3b, as initial conditions 
(see lNelson et al. 20101) . Figure [2{a) shows the temporal 
evolution of the hemispherical volume-average of toroidal 
magnetic field, [B l p\ 1 in the progenitor case over approxi- 
mately 5000 days, demonstrating the irregular cycles this 
model yields. Case S3 continues this behavior over about 
1300 simulated days starting from time t s . The tempo- 
ral evolution of [B^] in case S3 is shown in Figure EJb) , 
revealing two reversals of global magnetic polarity. 

Some caution should be used in interpreting any LES 
dynamo simulation, given the potential sensitivity of dy- 
namo action to magnetic dissipation and the nonlin- 
ear, nonlocal nature of turbulent magnetic induction, 
which makes reliable SGS modeling difficult. However, 
we believe the essential large-scale dynamics exhibited 
in this simulation are robust and are largely insensitive 



to the SGS model. Indeed convective dynamo simula- 
tions with differing prescriptions for SGS diffusion ex- 
hibit similar large-scale magnetic struc tures (B10; Bll; 
iGhizaru et al. 20101 : [Racine et al. 20"TTI ). 

Here wc will discuss buoyant magnetic structures which 
coherently rise above O.9Oi?0 while remaining connected 
to the large-scale toroidal wreaths. Using these crite- 
ria, we have identified nine buoyant magnetic loops, in- 
dicated by hash marks in Figure HJb) . Eight loops are 
seen in the northern hemisphere and one in the southern 
hemisphere. We expect that the apparent asymmetry is 
simply the result of having studied only two magnetic 
cycles. 

3. BUOYANT MAGNETIC LOOPS 

Buoyant magnetic loops arise from the cores of toroidal 
magnetic wreaths near the base of the simulated do- 
main. These wreaths have significant (B^) components 
that peak around 5 kG while also having strong non- 
axisymmctric fields. Figure HJb) shows a typical B^ 
configuration involving a negative polarity wreath in the 
northern hemisphere spanning 95° in longitude and a 
positive polarity wreath in the southern hemisphere ex- 
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tending over 270° in longitude. As demonstrated in cases 
D3 (BIO) and D5 (Bll), these magnetic wreaths arc 
highly nonuniform and display significant internal vari- 
ation as well as a high degree of connectivity with the 
rest of the domain. In case S3 portions of the wreaths 
can have coherent cores in which B^ can regularly ex- 
ceed 25 kG and have peak values as high as 54 kG. In 
these cores, bundles of magnetic field lines show very 
little local connectivity with the rest of the domain or 
even the other portions of the wreath. A single wreath 
of a given polarity may not form a coherent core at all 
or may have more than one core, and a single core may 
produce multiple buoyant loops. Of the nine buoyant 
loops investigated here to rise past 0.90i? Q , one coher- 
ent core produces four buoyant loops, another produces 
three, and two more cores each yield a single buoyant 
loop. 

Some of the coherent wreath cores can become buoy- 
ant magnetic loop progenitors or proto-loops. In these 
proto-loops the strong Lorentz forces result in highly sup- 
pressed convective motions. If we examine extended re- 
gions in the cores of wreaths with a local ratio of mag- 
netic to kinetic energy above a fiducial value of 100, we 
identify at least 35 proto-loops at the times where the 
nine buoyant loops arise. Thus the large majority of 
proto-loops do not evolve into mature buoyant loops, 
generally due to unfavorable interactions with convec- 
tive flows. When magnetic field strengths exceed 35 kG 
the proto-loops become significantly underdense as mag- 
netic pressure displaces fluid, causing buoyant accelera- 
tion. With some rise a proto-loop can enter a region of 
less suppressed giant cell convection. These flows will 
advect portions of the proto-loop downward at cell edges 
and upward in the core of the giant cells. The rise of the 
top of a magnetic loop is shown in cross-section by sam- 
pling Bj, roughly every 2 days in Figure Eta). Not all 
proto-loops become buoyant loops by our criteria. For 
example, loop B in Figures [2Ia,c) begins to rise but is 
prevented from passing 0.88 Rq when the top of the loop 
encounters a strong downflow. 

One way to track these buoyant loops is to use 3-D 
tracings of ma gnetic field lines u sing the VAPOR soft- 
ware package (jClvne et al. 20071) . In our simulations 
with finite resistivity, individual field lines do not main- 
tain their identity in time. However, one can achieve 
some measure of consistency as the structure moves and 
evolves by tracking field line ensembles. We track the 
very strong fields at the bottom of the loops near the 
base of the domain and we randomly seed large numbers 
of field lines (here 1000) in those highly magnetized foot- 
points at each time step. Figurc[3jb) shows a 3-D render- 
ing of magnetic fields lines for two sample loops near the 
peak of their rise. Similar field line tracings have been 
studied at various times during the rise of these loops. 

At maximum rise, the sample magnetic loop A ex- 
tends from 0.73-R© to O.93i? . The magnetic fields ex- 
ceed 40 kG at the base of the loop but become much 
weaker near the top of the loop, with field strengths 
as low as 2kG. Such loops are embedded in the much 
larger wreaths which have an average cross-sectional area 
of 13800 Mm 2 . The cross-sectional area of loop A is 
120 Mm 2 at 0.795i? Q and 520 Mm 2 at its peak radial po- 
sition of O.923i?0. Accounting for the continued expan- 



sion that would likely occur if this loop were able to rise 
further, the cross-sectional area is reasonable compared 
with the typical area of a large sunspot at the solar sur- 
face, which is roughly 2500 Mm 2 (jZwaan 19871) . If the 
loop were rising adiabatically over the same interval, the 
cross-sectional area should change in inverse proportion 
to the change in background pressure, which decreases 
here by a factor of 17.1, rather than the observed ex- 
pansion by a factor of 4.3. The top of the loop must 
then have a net outflow of heat or material in order to 
avoid expanding adiabatically. The loops show a mea- 
surable deficit in density and thermodynamic pressure 
relative to their surroundings, but they do not possess 
any detectable signature in temperature or entropy. This 
indicates that they are thermally "leaky" and able to 
equilibrate quickly compared to the timescale for radial 
motion. A simple estimate of the thermal diffusion time 
across one of these structures at mid-convection zone is 
on the order of 50 days, implying that there is likely 
also a divergent flow at the top of the loop, moving fluid 
along field lines. We see some evidence for such flows 
with roughly 1ms -1 speeds. 

Once a loop has begun to rise, its radial motion is 
dominated by advection and magnetic buoyancy. Fig- 
ure [UJe) illustrates the motion of loop A which begins 
to rise buoyantly at if,, while also indicating the com- 
ponents of the motion due to advection and magnetic 
buoyancy. To compare motion due to magnetic buoy- 
ancy, we define a magnetic buoyancy velocity w m b at 
the times sampled in Figure EJc). Magnetic buoyancy 
acceleration is here the fractional density deficit in the 
loop compared to the average density of the surround- 
ing fluid times the local gravitational acceleration. For a 
magnetic structure in local thermal equilibrium, this re- 
duces to the ratio of magnetic pressure inside the loop to 
thermodynamic pressure in the surrounding fluid times 
gravitational acceleration. To compute v m ^, we integrate 
the magnetic buoyancy acceleration over the intervals 
between times plotted in Figure (3{c) (roughly 2 days) , 
which likely provides a lower bound on this velocity. The 
advective velocity u ac j is the volume-averaged velocity of 
the surrounding fluid. The pressure and velocity of the 
surrounding fluid are calculated by taking averages over 
the convective updraft while excluding regions with field 
magnitude greater than 4 kG. Initially the sample proto- 
loop experiences an upward v m h = 46.1ms" 1 . After 3 
days of movement dominated by magnetic buoyancy, the 
loop gets caught in a convective updraft and v a d becomes 
greater than t> m b. Even though advective motions dom- 
inate, magnetic buoyancy continues to drive an average 
upward motion at 32.3 ms" 1 relative to the surrounding 
fluid. Continued buoyant acceleration of the loop as the 
magnetic pressure weakens is achieved because its den- 
sity perturbation decreases at roughly the same rate as 
does the background density stratification. Once the top 
of the loop has entered the main convective upflow it ex- 
periences advection at an average velocity of 53.1 ms" 1 . 
The presence of magnetic buoyancy forces allow this loop 
to rise in 14.6 days while the average upflow traverses the 
same distance in 21.7 days and magnetic buoyancy alone 
would require 30.6 days. 

Additional accelerations are present but not shown, 
including thermal buoyancy, which is significant early 
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in the rise of the loop, and magnetic tension, which is 
of the same order of magnitude as the advective mo- 
tion near maximum radial extent at 14.6 days and helps 
tether the loop to that height. Thermal buoyancy is dis- 
tinguished from magnetic buoyancy by averaging over 
the convective updraft but excluding regions with mag- 
netic fields above 4kG. An additional apparent motion 
at early times is produced as toroidal magnetic field used 
to track the loop is converted to radial magnetic field in 
the sides of the loop. Because advection plays a crucial 
role in the transport of these magnetic loops, their size 
scale is set by the size of the convective giant cells. The 
nine loops studied here have an average extent of 15.4° 
in longitude when measuring across the bottom of the 
loop, whereas the average distance between convective 
downflows in the equatorial region is 16.4° in longitude. 

4. REFLECTIONS 

In this paper we have presented a 3-D MHD simulation 
that combines turbulent convection, rotation, and strat- 
ification to produce solar-like differential rotation and 
wreaths of large-scale toroidal magnetic field at the base 
of the convection zone. These undergo cycles of magnetic 
activity and reversals of global magnetic polarity. Most 
notably the wreaths also exhibit buoyant magnetic loops 
capable of coherently traversing much of the convective 
layer. Such loops can only be realized when the field 
amplitude in a portion of a wreath exceeds 35 kG, the 
diffusion timcscale across the proto-loop (here 50 days) 
is much longer than the timescale for rise due to mag- 
netic buoyancy, and the interactions between rising loops 
and convective flows are favorable. These buoyant loops 
which appear at cycle maximum can have toroidal field 
strengths of 45 kG at their base and 5kG at their top. 
Their size scales are set by the size of the convective giant 
cells and they have cross-sectional areas at O.9Oi?0 that 
arc reasonable compared to the area of a large sunspot. 

We must be cautious in suggesting that these rising 
magnetic loops can make it through to the surface of 
the star. Our global simulations here only extend to 
0.97-R© and currently place an impenetrable boundary 
there, for we cannot cope with the intense small scales 
of convection seen as supergranulation and granulation 
near the surface. The presence of the domain boundary 
deflects all flows, leading to some uncertainty about the 
fate of the rising loops that could only be resolved by 
linking flows and magnetism in the upper reaches of ASH 
to another high-resolution compressible domain closer to 
the surface. This is a task we are now pursuing in parallel 
with global modeling. 

It is noteworthy that within this simulation convection 
generates differential rotation which in turn generates 
toroidal flux which then buoyantly destabilizes and rises. 
Each link in this chain is physically well established. 
Our primary accomplishment here is to capture all these 
processes self-consistently within a single simulation. 
This represents an essential step toward unifying nu- 
merical models of global-scale convective dynamos and 
surface flux emergence. 
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